SAT in P 


Abstract. There are a large number of papers that claim that there are problems that once 
solved lead to an efficient solution of a wide range of problems, classified as NP. In this 
paper we will not only question the existence of this class of NP-co problems, but we will 
also explain their limitations in engineering and give a polynomial-time solution to SAT, 
one of these emblematic problems. The resolution will be so trivial that it will even be 


possible to practice it on paper. 


1. Introduction 


A Turing machine is a mathematical model of computation that defines an abstract 
machine that manipulates symbols on a strip of tape according to a table of rules. lt is 
considered that in that machine it is possible to denote any enumerable problem, and 
when each of the intermediate states are exactly determined Alan Turing designated that 
as an automatic Turing Machine [9] (a-TM, deterministic TM), if, on the contrary, it is 
always possible that the next state can be a finite number of possible states then that 
machine is denoted a choice Turing Machine (c-TM, non-deterministic TM, NDTM). 


The Boolean satisfiability problem (sometimes called propositional satisfiability problem 
and abbreviated SATISFIABILITY or SAT) is the problem of determining if there exists an 
interpretation that satisfies a given Boolean formula. That problem is usually reduced to 
formulas which are a Boolean product of sum of literals. Where a literal is a Boolean 


variable negated or not, and a sum of literals is called a clause. 


If we write on a tape the expression of a product of clauses and ask what values the 
variables must take to satisfy the formula, if we had those values we would only need 
polynomial time to determine if the formula is satisfied. That is why in a NDTM it is 
considered to be bounded in polynomial time, because in the absence of determining a 


finite set of values proportional to the input it would be solvable in polynomial time. 


Therefore, the SAT problem is considered trivially within the NP class. Having a 
mechanism that allows us to know what values the variables should take in polynomial 


time would convert SAT to the class P. 


Method used 


Para aseverar la afirmación usaremos dos métodos: uno sencillo de comprender y otro 
que exige un cierto compromiso matemático. El segundo método tiene la ventaja de que 
funciona más rápido cuantas más máquinas estén implicadas debido a que es fácil de 


distribuir en una red de computadoras que trabajen en paralelo. 


Method 1 


Para el primer método dispondremos de cualquier fórmula bien formada y encontraremos 
la manera de expresarla como un sistema de ecuaciones lineal módulo 3 con una única 
restricción: que todas las variables deberán ser distinto de cero. 

Es fácil encontrar un sistema de inecuaciones módulo P, siendo un número primo impar, 
que represente un producto de sumas de literales booleanos. Por eso, como la restricción 
de que las variables no sean O equivale a encontrar un equivalente en formato de 


inecuación, se ha elegido la representación más sencilla. 


Empecemos por estudiar una ecuación sencilla: x + y + A+ B =1 (mod 3) sabiendo que x, 


y, A, B no puede ser 0. 























x+ y + A+ B= 

1 1 1 1 1 
1 2 2 2 1 
2 1 2 2 1 
2 2 2 1 1 
2 2 1 2 1 




















Nótese que son cinco entradas, si consideramos que x mod 2 es la manera que tenemos 
de establecer la equivalencia hacia la aritmética modular, entonces 2x = 2x. Además 
observamos que B = xy, salvo en la última entrada; lo cual es un operador interesante de 
representar. Más en concreto, B =x O y (mod 2). Así que si estudiamos qué pasa con los 


negados de x e y entonces ocurrirá que B ==x O 7y y creamos la siguiente tabla: 
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2x + 2y + C+ B= 
2 2 2 1 1 
2 1 2 2 1 
1 2 2 2 1 
1 1 1 1 1 
1 1 0 2 1 




















Observamos ahora que cuando combinamos las ecuaciones 
x+y+A+B=1 (mod 3) 
2x + 2y + C+B=1 (mod 3) 
conseguimos hacer que B=xy equivalga a una operación XOR, y además vemos que A 
corresponde con la operación AND y C corresponde con la NOR. Por tanto ya tenemos 
operadores más que de sobra para representar cualquier ecuación booleana mediante un 


sistema de ecuaciones en base 3, con la restricción de que sus variables no se anulen. 


Una de las operaciones que se debe tomar en cuenta a la hora de resolver nuestro 
sistema de ecuaciones es cuándo la solución ya está determinada, y no es posible 


encontrar grados de libertad. Las posibilidades válidas que se tienen son sólo dos para 
una variable: 


y otras cuatro para dos variables: 
x + y =1 (mod 3) <=> x =2 (mod 3) £ y =2 (mod 3) 
x + y =2 (mod 3) <=> x=1 (mod 3) 8 y =1 (mod 3) 
x + 2y = 1 (mod 3) <=> x =2 (mod 3) 8£ y =1 (mod 3) 
x + 2y =2 (mod 3) <=> x=1 (mod 3) 8 y =2 (mod 3) 
También existen 4 ecuaciones homogéneas simplificables: 
x+y+z=0 (mod 3) <=> x + 2y =0 (mod 3) 8 y + 2z =0 (mod 3) 
x + 2y + Z=0 (mod 3) <=> x + y =0 (mod 3) 8 y + z =0 (mod 3) 
x + y + 2Z =0 (mod 3) <=> x + 2y =0 (mod 3) 8 y + z =0 (mod 3) 
x + 2y + 2zZ =0 (mod 3) <=> x + y =0 (mod 3) 8 2y + z=0 (mod 3) 
Sabiendo que, dependiendo de cómo se quiera implementar el solver, al encontrar 


ecuaciones como x + By = 0 (mod 3) automáticamente se pueda eliminar la columna y, así 
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como esta ecuación para que, una vez resuelto el valor de x se deduzca: y = 2Bx mod 3. 
A las aclaraciones hechas sobre simplificación y determinación de ecuaciones las 
llamaremos clausura del sistema: en el sentido de que es posible que haya que 
aplicarlas reiteradamente hasta que encontremos un sistema de ecuaciones con 
suficientes variables por ecuación como para continuar resolviendo de manera 


determinística. 


Para demostrar que es posible crear un solver que resuelva en tiempo y espacio 
polinomiales el problema de SAT empezaremos por recomendar que para cada par de 
literales a los que debamos hacerle una operación lógica se le calcule las cuatro 
ecuaciones para las variables ¡, j coni<j: 
V¡+ V; + Ars + X¡ = 1 (mod 3) 

V¡+ 2Vj¡ + Aro + 2X¡ = 1 (mod 3) 

2V¡+ V; + Aos¡¡ + 2X¡ = 1 (mod 3) 

2V¡+ 2V; + Apo + X¡ = 1 (mod 3) 
Nótese que a la hora de formar la matriz las columnas que representan a V, estarán bajo 
un orden, las de A en otro y las X en otro, así que dispondremos las columnas de índice 
más pequeño con las X, después colocaremos las A y, finalmente pondremos las V en ese 
orden. 
Además, se observa una condición adicional: 

Ar + Aro + Aor + Avo¡ = 1 (mod 3) 

A esta submatriz de cinco nuevas filas y cinco nuevas columnas a añadir por cada 


relación ¡, j que se encuentre lo llamaremos la caja i, j. 


De esa manera, una vez formada la matriz, si se busca una base a partir de las primeras 
columnas lo más probable es que se dejen las variables V, como parámetros según 
corresponda. 

Otro aspecto a tomar en cuenta es qué pasaría si la matriz acaba teniendo más filas que 
columnas, o la misma cantidad; en ese caso ya se tendría la solución. Y si alguna de las 
variables tuviera que ser O o acaba en una inconsistencia eso quiere decir que no hay 
solución. 

Por supuesto, como ocurre con la triangulación de Gauss, si se descubre que una fila 


acaba en O = 0 entonces deberá ser eliminada. Así que procedemos a mostrar el 


algoritmo: 






































1. Se transforma la fórmula booleana en un circuito O(n?) 
2. Se ordenan las columnas que serán necesarias para el circuito O(n) 
3. Se genera una caja por cada puerta del circuito O(n) 
4. Se encuentra la identidad de la matriz con las primeras columnas O(n?) 
5. Se calcula la clausura del sistema O(n) O(n?) 
6. Si no hay más columnas que filas: Se resuelve y FIN O(n) 
7. Se lanza un supuesto desde las últimas columnas O(n) 
8. Si el supuesto acierta: Se resuelve y FIN O(n) 
9. Si el supuesto falla: se añaden las cajas que niega el supuesto O(n) 
10. Se reducen las nuevas filas a partir de la identidad (normalizar)  |O(n) 
11. Si las nuevas filas fueron eliminadas: Insatisfacible y FIN. 1 

12. Salta al 5 O(n) 

















Nótese que se salta a 5 desde 12 tantas veces hasta que el número de filas sea igual o 
superior al número de columnas, sabiendo que una caja añade siempre cinco filas nuevas 
para cinco columnas, salvo que el par del circuito ¡,| ya haya sido usado antes, en cuyo 
caso se incorpora sólo la caja del circuito que no ha sido referenciada. 

Esa caja anulará los parámetros sobre la base anterior para hacer que la próxima vez el 
supuesto se haga con menos grados de libertad. Por lo que por cada supuesto se reduce 


en al menos una el número de columnas que generaron la indeterminación. 














Base A V 
0 0 Xx Xx Xx x 
0 1 0 Xx Xx Xx x 
0 0 1 Xx Xx Xx x 




























































































Base A V 
1 0 0 Xx Xx Xx x 
0 1 0 x Xx Xx x 
0 0 1 x Xx Xx x 
0 0 0 C C 0 0 
0 0 0 C C 0 0 
Base A V 
1 0 0 0 0 Xx x 
0 1 0 0 0 Xx x 
0 0 1 0 0 Xx x 
0 0 0 1 0 0 0 
0 0 0 0 1 0 0 





























Entre las columnas nuevas, las filas nuevas, las columnas que se aprovechan y que 
forman parte de la nueva base, al final se van arrinconando los parámetros libres y se va 
a aumentando el número de filas a satisfacer. Pero nunca podrá aumentar más allá de las 


columnas que vaya incorporando a la base. Por eso la iteración es lineal. 


Method 2 


It is possible to speak of a third method whose complexity is lower, provided that it takes 
advantage of having several machines running in parallel. Using matrix notation, n 
machines on an input of n clauses can distribute the work of multiplying the matrix by itself 


and go from having a cost of O(n? log n) to being O(n log n) for large inputs. 


To do this, the following rules must be available: 

1. Two clauses of size n and m generate a submatrix of 1's of nxm. 

2. Every literal ¡ and j that are opposites becomes 0 in (i, j). 

3. Every literal ¡ and j that are the same becomes 0 in k 41 (i, k) and k 4j (k, j). 


In this way a matrix is formed with the submatrices relating all clauses to themselves to 
form a square matrix of size mxm where m ¡is the sum of all literals along the input, or the 


product of the number of literals per clause times the number of clauses if it was the same. 


Partiendo de estos elementos definiremos un cuerpo mediante los operadores de 


producto y suma de la siguiente manera: 


AxB = A8((AB>=n) | |) 
where n is the number of clauses and | identity matrix. 
A+B =AxorB xor | 


elemento neutro | 


Para proceder a demostrar la propiedad distributiva nos valdremos de lo que se 
demuestra en los apéndices, sabiendo que como el operador (A(B € C)>=n) se puede 


traducir como: 


(XA¡ 8(Bx € Cy)>=n) entonces se deduce trivialmente que equivale a 
(Aj 4Bx O Aj¡8Ci)>=n) para garantizar así la propiedad distributiva entre la operación 


matricial y el xor. 


The matrix product consists of calculating the scalar product of row and column and 
returning 1 if the result is greater than or equal to n; n is the number of clauses. 

When the matrix multiplied in this way with itself gives the same matrix we will say that it 
has reached a Noetherian ideal, and if it maintains the main diagonal then we will proceed 
to find a k-th solution. To find the k-th solution we only have to choose a pair of clauses 
whose submatrix has more than one 1 and we select only one 1, then we calculate ¡ts 
noetherian ideal and if the result is O then it was inconsistent, if on the contrary we repeat 
the procedure until the matrix reflects a single minterm. 


For the next fórmula 




































































TE (x Yo ETA A REA EZ) (ex +y +2) 
x 1 0 0 0 1 1 0 0 1 
+ y 0 1 0 1 0 1 0 1 0 
+ 2) 0 0 1 1 1 0 1 0 0 
(ax 0 1 1 1 0 0 1 0 0 
+ =y 1 0 1 0 1 0 0 0 0 
+ z) 1 1 0 0 0 1 0 0 1 
(ax 0 0 1 1 0 0 1 0 0 
+ y 0 1 0 0 0 0 0 1 0 
+2) 1 0 0 0 0 1 0 0 1 
TONTA STATES 
0 1 1 10:00: 1. 1 0 1 1|/1 O 0| [O O 1|/1 0 0 
1.0 1AO 1 O[[1 O 1|+|1 O 1/0 1 O[+[O 1 OJO O O|[>3 
1.1.0 0 0 1/1 1 0 11 0/0 O 1 1.0 0/10 O 1 
0 1 1 0 1 1 0 1 1 0 0 1 
1.0 1jAp[1 O 1|+[1 O 1|+[0 O 0|[>3 
1.1.0 1.1.0 1.10 1.0.0 


0 1 1 0 0 1 0 0 1 
1 0 110 0 OlÍ=I|0 O O|=xA7ZVxXAZ 
11.0 10.0 100 


What this algoritham should induce us is that it is still possible to find methods that on the 
most powerful machines can give us more efficient results than the quadratic one. 
It is easy to see this algorithm will work because this little lemma: 

4/(4,-1)=0>] | 4,=u(-n+) 4,) , where u(x)=1 iff xis not negative, else 0. 
Considering a matrix is formed by subtables, each table with only one solution is enough 
to be a sufficient condition to ensure the satisfiability of the complete formula under a 
Noetherian ideal because every solution will get catched and: 














SAT 
SAT(f) > Noetherian(f) 
1 1 1 
0 1 1 
0 1 0 

















So after finding a solution if there is noone f is unsatisfiable. 
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Apéndices 


Lema 1. 
















































































Lema 2. 











































































































Lema 3. 
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Lema 4. Propiedad Distributiva adulterada de O y V. 
(XOY)VZ=(XVZOYVZ)VZ 

Dem. 

Caso Z = 1: 1=1 

CasoZ=0:(X6Y)VO=(XVO0OSYVO)VO 


Demostrado. 


Si para formar un cuerpo tenemos un operador intermedio $ que cumple : 
1) existe un elemento e donde AGVe=e8A=A 
2) se cumple quee VA=AVe=A 
3) el producto A x B tiene la forma AA(A OS B) V e, 
4) la suma A + Btiene la forma AB Oe 
5) y se cumple A9G(BOC)=ASBO ASC, propiedad distributiva con el O 
entonces se cumple: 
Ax(B+C)=AxB+AxGC, la propiedad distributiva 
siempre que A, B y C tengan un factor común 


Dem. 


En el proceso se va a cambiar la expresión a una expresión siempre equivalente 


1. AxB+AxC 


12 


2. (AN(A0B)Ve) S(AN (AO C)Ve) Oe 
3. (AA(AOB)ODAN (A9C)) Ve 

4. AN(A0B)O (AQ C)) Ve 

5. AN(A9 (BOC) Ve 

6. AMAO(BOC)IO e) Ve 

7. ANA0((B+C)0 e)0 e)Ve 

8. (AVeJA((A9 ((B+C)0 e)O e) Ve) 

9. AM((Z,0 e)Ve) Z,=A09((B+C)O e) 
10. AAM((Z¡Ve0 eVe0de) 

31.AM((AO ((B+C)O e)) Ve) 

12. AA((A0Q(B+C)O A)Ve) 
13.AM((A0(B+C)Ve)O (AVe)O e) 

14. AM((A0Q(B+C)Ve) GANAVe)O AñNe 


15.Ax (B+C)+A 


[Por definición de producto y sumas] 
[Por lema 1 aplicado en paso 2] 
[Por lema 2 aplicado en paso 3] 
[Por propiedad distributiva] 

[Por lema 3 aplicado en paso 5] 
[Definición +] 

[Distributiva A, V] 

[Propiedad neutra e con V] 


[Lema 1 aplicado en paso 9] 


[Lema 2 sobre 11] 
[Lema 1 sobre 12] 


[Lema 2 sobre 13] 


Supongamos que queremos estudiar qué pasaría si los elementos fueran múltiplos de un 


factor común. 
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AxDx(AxB + AxC) 


AxDx(Ax(B+C) + A) 


AxDxAx(B+C) + AxDxA 


AxDxAxB + AxDxAxC +AxDxA + AxDxA 


AxDxAxB + AxDxAxC 


